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ABSTRACT 


An empirical orthogonal function (EOF) analysis is used 
to represent environmental wind fields associated with 
selected western North Pacific tropical cyclones from 
1979-1983. Composite synoptic patterns within five  past- 
motion categories of cyclones are studied using 5, 10, 20 
and 35 EOF modes. Significant differences in the composite 
fields are found between categories, which suggests that the 
wind EOFs are capable of representing synoptic patterns 
associated with tropical cyclone motion. Regression equa- 
tions are developed to post-process the forecast track of 
the One-way Tropical Cyclone Model (OTCM). The predictors 
for the regression analysis include the EOF coefficients, 
past motion, Storm position, date, intensity, backward 
extrapolation and forward track displacements. For the 
dependent sample, the mean /2-h track forecast error for the 
modified OTCM is 383 km, which represents a 210 km improve- 
ment relative to the unmodified OTCM. Thus, the regression 


scheme has great potential for operational application. 
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AUN TRODUGTEON 


A. BACKGROUND 


One of the most difficult problems tier op lca) 
meteorology is to forecast the movement of tropical 
cyclones. Although operational tropical cyclone track 


forecasting remains rather subjective, a number of objective 
forecast aids provide guidance as to the future storm track. 
No one forecast aid has yet provided consistent guidance in 
the wide variety of observed weather situations. Each aid 
has strengths and weaknesses. Better tropical cyclone 
forecasts hinge on a better utilization of existing forecast 
aids (Neumann and Pelissier, 1981; Tsui, 1984). 

In this study, synoptic information is used to develop a 
Statistical-synoptic forecast technique which will be used 
` to modify an existing forecast aid used by the Joint Typhoon 
Warning Center (JTWC) in Guam. The modification scheme (see 
Chapter 4) should enable the JTWC forecaster to better 
utilize the existing aid. It is appropriate to give a brief 


overview of the forecast aids that use synoptic information 


to make cyclone track forecasts. The original 
statistical-synoptic techniques (e.g., Riehl et al., 1956) 
used latitudinal and longitudinal differences of 


geopotential values to estimate the steering components of 
the storm. Two examples in use today by the National 
Hurricane Center (NHC) are the NHC67 and NHC72 models. 
These models use the current and 24-h old 1000, 700 and 500 
mb geopotential height data to modify a preliminary forecast 
based on climatology and presistence (Neumann and Pelissier, 
T981). In these models, the most significant geopotential 
heights, which are positioned across the storm, act as 
steering predictors. 

Statistical-synoptic techniques typically use a grid 


that is relocated each forecast period so that the storm is 
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always at the same gridpoint. Synoptic data, suchi: 
geopotential height values, are subjectively or objectively 
interpolated onto the grid. The gridpoint values are then 
used to develop regression equations to predict zonal and 
meridional displacements of the storm. The regression 
analysis typically chooses only a few key gridpoint values. 
Therefore, this approach incorporates only a small portion 
of the original synoptic field. 

The synoptic-scale data used in statistical-synoptic 
models have not been limited only to analyses at the 
forecast time. The NHC73 model uses numerically predicted 
values (from a global model) as potential predictors 
(Neumann and Lawrence, 1975). The NHC73 model has proven to 
be one of the most consistent methods in the Atlantic region 
(Neumann and Pelissier, 1981). 

Synoptic-scale information is also used in some 
. analogue-type prediction schemes. Analogue techniques are 
based on the assumption that a tropical cyclone will move 
similar to previous storms that occurred within some spatial 
region and temporal interval of the current storm. A 
computer algorithm searches numerous cases until a family of 
storms is found with the same general characteristics as the 
current storm. The storm is then forecast to move along the 
weighted mean track of the family of storms. An example of 
an analogue routine that uses synoptic information was 
described by Jarrell and Somervell (1970). They included 


the location of the subtropical ridge and midlatitude trough 


in their routine. A weakness in these analogue techniques 
is the failure to give a forecast under anomalous 
conditions. 


In all of the statistical techniques discussed above, 
the synoptic field has been represented by gridpoint values. 
Only selected points from the field, or a characterization 


of the field such as the latitude of the trough line, have 


been used. Adjacent gridpoints in the field may be highly 
correlated, and thus not provide independent information. 

An alternate approach to a gridpoint representation of a 
Synoptic field has been described by Shaffer (1982). He 
used empirical orthogonal function (EOF ) analysis to 
represent the 500 mb geopotential height fields on a grid 
centered on the tropical cyclone. Shaffer and Elsberry 
(1982) demonstrated that the coefficients from the EOF 
analysis could be used as synoptic forcing predictors in a 
statistical-synoptic track prediction scheme. 

In a similar study, Wilson (1984) used EOF analysis to 
represent the 700, 400 and 250 mb wind component fields on a 
grid centered on the tropical cyclone. Wind fields may have 
greater utility than geopotential height fields in the 
tropics because the geopotential height gradients are weak 
and the flow is far from geostrophic. This is especially 
. true for smaller scale motions such as tropical cyclones. 
Wilson demonstrated that the coefficients from the wind EOF 
analysis could be used as synoptic forcing predictors in a 
statistical track prediction scheme. 

The numerical model is another type of track forecast 
Scheme that uses synoptic information. One of the first 
numerical models was the SANBAR barotropic model, which was 
based on the assumption that momentum advection is the 
primary physical mechanism for motion of intense tropical 
cyclones (Sanders and Burpee, 1968). This simple 
single-level model steers the storm based on the 
Aeee Cale i vertically integrated current in which it is 
embedded. 

Recently, several meteorological agencies have 
successfully applied fine-resolution or nested numerical 
models to the problem of predicting tropical cyclone 
movement. The National Meteorological Center (NMC) uses a 
Moveable Fine Mesh (MFM) model. The fine mesh grid is 


relocated relative to the storm center during a forecast. 
The MFM has 10 vertical layers and includes physical 


processes Similar to the primitive equation models used in 


mid-latitudes. The U. T Navy's Fleet Numerical 
Oceanography Center (FNOC) has two numerical models for 
predicting tropical cyclone movement. The One-way Tropical 


Cyclone Model (OTCM) is discussed in detail in Chapter 4. 
Although the OTCM is a three-layer, primitive equation model 
(Hodur and Burk, 1978), it has proven to be one of the best 
forecast aids used by the JTWC (Tsui, 1984; Peak and 
Elsberry, 1986). The other dynamical tropical cyclone model 
used by the Navy is the two-way interactive Nested Tropical 
Cyclone Model (NTCM). The model uses a fine resolution grid 
(41 km) near the tropical cyclone and a coarse grid (205 km) 
away from the cyclone. The initial results of the NTCM were 
described by Harrison (1981). 

Numerical weather prediction in the tropics has not been 
as successful as mid-latitude numerical prediction for a 
number of reasons. Tropical regions have very poor data 
coverage compared to the mid-latitude regions. This results 
in a poorly defined initial field for the tropical model. 
Another problem is that the physical processes involved in 
tropical circulations are not fully understood. | 

In summary, Several forecast schemes use synoptic 
information to forecast the movement of tropical cyclones. 
These forecast models range from statistical to numerically 
based schemes. The major weakness of statistical schemes is 
forecasting for anomalous conditions, while dynamical 
forecasting tends to exhibit systematic errors. Since the 
systematic errors are repetitive, a statistical 
post-processing scheme (see Chapter 4) may be used to reduce 


these errorss 
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DOEUBJECTIVES 

The current tropical cyclone forecast requirements for 
JTWC are summarized by  Sandgathe (1985). Accurate 48-hour 
forecasts are normally required for air and naval base 
evacuation and preparation. Accurate 72-hour forecasts are 
required for ship routing and major operational exercises. 
Better accuracy in the forecasts is desirable to avoid 
costly multiple evacuations of aicraft or overloading of 
strategic bases. Guided by the current requirements at 
E this study will eoncentrate on improving the accuracy 
of forecasting the cyclone position “at the 72-h interval. 
One of the primary forecast aids at JTWC is the One-way 
Tropical Cyclone Model (OTCM). A method that can help the 
JTWC forecaster better utilize the OTCM would be beneficial. 
Therefore, a goal of this Study is develop a modification 
scheme for the OTCM. Any improvement in the OTCM forecasts 
Should be beneficial, since it has proven to be one of the 
best forecast aids used by JTWC (Tsui, 1984; Peak and 
]sberry. 1986) 

The same data set used by Wilson (1984) is used in this 
study. It consists of EOF coefficients based on synoptic 
wnd information, along with the appropriate tropical 
cyclone information. A detailed description of the data set 
is found in Chapter 2. The first step in this research is 
to demonstrate that the wind-based EOF coefficients can be 
used to define synoptic effects acting on tropical cyclones. 
One advantage of EOF analysis (see Chapter 2) is a 
considerable reduction in the storage required to represent 
the synoptic forcing, as the entire wind field surrounding 
the tropical cyclone can be defined by a small set of EOF 
coefficients. A compositing technique using the EOFs to 
define the synoptic forcing of the tropical cyclone track is 
presented in Chapter 3. The second step in this study is to 


utilize the synoptic description in terms of the wind EOFs 
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to adjust the OTCM (see Chapter TE E AwSsUccessimi 
post-processing scheme should provide the JTWC  forecaster 
with guidance as to when and how the OTCM forecast should be 
modified. Finally, Chapter 5 addresses the implementation 
of the OTCM post-processing technique and suggests future 
research related to this study. 

To Summarize, there are two objectives for this study. 
First, demonstrate that the EOF coefficients can be used to 
differentiate between synoptic Situations affecting tropical 
cyclone motion. The second objective is to derive a 


Statistical post-processing scheme for modifying the OTCOM. 
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AIDA TA ACQUESTIION AND"FIELD DEFINITION 

This study employs the same data set used by Wilson 
(1984). Western North Pacific tropical cyclone information 
from 1979-1983 is obtained from the annual tropical cyclone 
reports published by the Joint Typhoon Warning Center (JTWC) 
in Guam. The data set includes warning and best-track 
positions at six-hour intervals, as well as the estimated 
intensity of the storm (maximum sustained winds). 

The U. S. Navy's Fleet Numerical ELTE Center 
provided operational wind fields from the Global Band 
Analyses (GBA). The GBA wind fields are produced every 12 h 
on a Mercator projection true 2c? 7.5" NE The grid 
resolution is 2.5? longitude by approximately 2.5? latitude. 
The GBA provide complete longitudinal coverage between 41?S 
and 59.8°N. When a tropical cyclone is present, eight bogus 
winds at the surface are inserted at 80 km from the center 
of the cyclone, and vertical coupling is achieved via the 
thermal wind relation using temperature analyses at 
intermediate levels. Upper-air observations include 
rawinsondes, pibals, aircraft and cloud motion vectors. T 
regions with no observations, the analysis will be a blend 
of the 12 h old analysis and 5% climatology. The interested 
reader can refer to the U.S. Naval Weather Service (1975) 
for a detailed description of the GBA. The GBA fields are 
available for the zonal and meridional wind components at 00 
andi l2 GMT since 1975. For this study, only data between 
1979-1983 at 700, 400 and 250 mb are used. 

Several criteria are imposed in selecting a storm for 
inclusion in the data set. A tropical cyclone must have 
matured to at least tropical storm intensity (maximum 


sustained winds of 18 m/s or greater) and must have been 
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aya 
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ae 2.l.. The relocatable RU Ee LE grid with a fixed zonal 
28 Que [orga Separation of 2/7.8 km The vortex center 


is located at gridpoint (16,9). from the lower left 
corner, and is indicated by the large dot. 


west of the dateline. The GBA must have been available for 
the zonal and meridional wind components at 700, 400 and 250 
mb. Finally, the JTWC warning positio AC OO cL elt 
must have been at a latitude south of 34.69?N. This larter 
restriction is necessary to ensure that the GBA winds are 
available for a sufficient latitudinal extent north of the 
cyclone center. A total of 1357 cases meet these criteria. 
The grid in this study is the same™as "Wilson s (1994) 
relocatable grid (Fig. 2.1). It consists of 527 data points 
with a fixed zonal and meridional separation of 277.8 km 
(Omen mie There are 31 zona lend 317 meridiora 


BRINDO: For each case, the grid s posibIOuede cibi 
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the tropical cyclone center is always located at the grid 
Paine (16,9). The warning position from JTWC is used to 
indicate the center of the tropical cyclone. For each case 
that meets the above criteria, the zonal and meridional wind 
components at 700, 400 and 250 mb are interpolated onto this 
equidistant grid from the GBA using a bilinear interpolation 
method (Wilson, 1984). 

RKecent ly, statistical stúdies On tropical cyclone motion 
have shown an improvement in 24-h forecast errors using a 
grid system oriented with respect to the initial storm 
heading (Shapiro and Neumann, 1984). However, the 
advantages for such a grid system are substantially 
diminished at the 48 to 72-h forecast times. Por enis 
reason, this study uses a geographically-oriented grid. 
However, five past-motion categories will be used to group 


similar storm tracks based on the prior 12-h motion. 


TE Bae EOR METHOD 

The method of principal component or empirical 
orthogonal function (EOF) analysis was first applied to 
atmospheric sciences by Lorenz (1956). In the past, the EOF 
analysis has been used as a map typing tool or as a 
technique for reducing the dimensionality and explaining the 
variance structure of a field (Wilson, 1984). For example, 
Stidd (1967) used EOF analysis to describe the variation in 
average monthly rainfall in Nevada. Brown (1981) divided a 
large sample into smaller discrete subsets by map typing 
based on the coefficients derived from the EOF analysis. 
Other examples include studies by Kutzbach (1967), Hardy and 
Walton (1978), Legler (1983) and many others. 

A one-dimensional data field can be generally 
represented by simple orthogonal functions such as sines and 
cosines. Similarly, the eigenvectors derived from EOF 
analysis can be used to reconstruct a two-dimensional data 


Treid. A major advantage in the EOF analysis is that the 


In 


eigenvectors are ranked in decreasing order according to the 
percent of the variance accounted tor (Morrison, 1967). For 
example, the first eigenvector explains the largest variance 
in the field. The variance unexplained by the first 
eigenvector would be the residual. The second eigenvector 
would then explain the maximum variation remaining in the 
residual field and so on. Since the first few eigenvectors 
describe a large portion of the total variance in a sample, 
one may approximate the total field by retaining only these 
few eigenvectors. As a result, the dimensionality of the 
data set can be reduced considerably, which makes the EOF 
method a cost-effective means of representing large data 
Sees. 

The procedures for the EOF analysis used in this study 
follow the treatment outlined by Wilson (1984). A scalar 
EOF analysis is performed using the zonal (u) and meridional 
(v) components at 700, 400 250 mb and a sample size of 682 
cases. From the 1357 cases that meet the desired criteria, 
alternate cases within a storm are chosen to make up the 682 
cases. At each pressure level and for each wind component, 
a 527 X 682 matrix A is formed such that the interpolated 
wind field for a case would appear as a column. The A 
matrix is normalized by subtracting the sample mean at each 
gridpoint, and then dividing by the sample standard 
deviation at that grid point. The normalized matrix Z is 


then defined with elements z(i,j) given by 


203) = LaS) DEA Sear (2.1) 


where a(i,j) are the elements in matrix A and b(i) and s(i) 
are the mean and standard deviation of the elements in row i 
of matrix A. In other words, b(i) and s(i) are the mean and 
standard deviation of the wind component for all the cases 


at a particular gridpoinom is The normalization of the wind 
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components in this way ensures that the  extratropical 
pegdens ofthe grid, where the variability is generally 
large, will not dominate the less variable tropical regions 
in the EOF analysis. A disadvantage in uSing a standardized 
matrix is that a small amount of smoothing may occur in the 
resultant eigenvectors (Kutzbach, 1967). 


The symmetric correlation matrix (R) can then be defined 


by 


m= ZZ /n; 2) 


where Z' is the transpose of Z and n is the number of cases. 


À scalar y can be defined by: 


pee Re, (2.3) 


where y is the correlation between an arbitrary vector e and 
the normalized data matrix Z. The vector e has the 
constraint that it must be normalized to length one (i.e., 
e'e = 1). It was shown by Morrison (1967) that (2.3) can be 


tug Pen in matrix form as 


EY = RE. (2.4) 


In this study, E is a 527 X 527 matrix, and the elements of 


Y are eigenvalues found by solving |R - YI| = O (where I is 

the identity matrix). Each column in E is an eigenvector e 

associated with a Single eigenvalue y. . D ict he 
i 


first eigenvector explains 


m 
ey 
Je 


of the total variance. The second eigenvector explains the 


largest amount of variance unexplained by the first and so 


on. 
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The normalized data matrix Z may be reconstructed by 


calculating the EOF coefficients or multiplicrsome( Stadd, 


1967; and others). These coefficients are found by solving 
C- E'Z, (259 
where C is am by n matrix (527 X 682 for this study)" The 


first column of C is the orthogonal coefficient vector 
corresponding to the first case, and so on. The normalized 


matrix Z can now be retrieved by 
Z= EC (22806) 


The jth case, stored in the jth column ot the Z matrix edm 


be represented by 


m 
z(1,3) s 9 eG kc S NETTEN Soy ( 2:898) 
k21 
where z(i,j) are the normalized gridpoints values for the 
jth case, e(i,k) are the elements of eigenvector k, c(k,j) 
are ene coefficients for the jth case and m is the number of 
modes (eigenvectors) retained. 


Since the elgenvectors in E are ranked according to the 


amount of variance they represent, a large portion of a 
particular field signal can be represented by a few 
eigenvectors. Wilson (1984) used the Monte Carlo approach 


(Preisendorfer and Barnett, 1977) to determine the number of 
elgenvectors needed to represent the signal of a field and 
eliminate unnecessary noise. The number of EOF modes that 
represents signal varies from 19 to 35 depending on the 
pressure level and wind component field (his Table VI). 
Wilson demonstrated that retaining only 35 of the 527 EOF 
modes (eigenvectors) explained between 81.7 to 93.2 percent 


of the uor v signal at the various pressure levels (his 
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Table VI). Mieretore only the first 35 EOF modes are used 
I this- study: 
Wilson pointed out that the eigenvectors computed with 


the dependent sample of 682 cases could be used with 


independent data. However, when a new case is added, the 
eigenvectors no longer exactly represent the maximum 
variation of the new plus the old observations. Wilson 


demonstrated that this does not introduce a significant 
error as long as the original correlation matrix R is large. 
The procedure for adding independent data is very similar to 
that outlined above. The A matrix is formed using the new 
cases. Each element in A is normalized by the gridpoint 
mean and standard deviation found from the original cases. 
The EOF coefficients are then calculated using (2.5), where 
the E matrix contains the eigenvectors found from the 
original cases. This procedure is used to determine the EOF 
. coefficients for the remaining 675 cases not used in the 


original EOF analysis. 


C. CT/AT COMPONENTS 

Numerous studies have been performed to determine the 
accuracy of tropical cyclone track forecasts (e.g., Jarrell 
et al., 1978; Neumann and Pelissier, 1981; Tsui, 1984; and 
Elsberry and Peak, 1986). Tsui studied the accuracy of 
various forecast aids used by the JTWC. He considered 
famine. track, speed and forecast errors. The forecast 
error was defined as the great circle distance between the 
forecast and the verifying position. The latter is based on 
the subjectively smoothed best track, which is produced in a 
post-storm analysis. reel error is defined as the 


right-angle distance from the forecast position to the best 


track. The other error definitions can be found in Tsui 
(1984). Neumann and Pelissier (1981) introduced another 
error measurement which involved cross-track (CT) and 
along-track (AT) components. Elsberry and Peak (1986) used 
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a similar method to calculate CT/AT components relative to 
an extrapolated track based on the present and past 12-h 
warning positions. 

In this study, CT and AT components are calculated with 
respect to a CLImatology and PERsistence (CLIPER) forecast 
track (Fig. 2 72E A western North Pacific version (called 
WPCLPR), provided by C. Neumann of the National Hurricane 
Center, is used (Xu and Neumann, 1985). The CLIPER track 
forecast is considered as a no-skill forecast, since the 


scheme does not use any synoptic information (Neumann and 


Pelissier, 1981). For each case, the CLIPER forecast is 
made using the -12 and -24 h warning positions, intensity 
and date of the storm. Warning positions are used because 


the best-track positions are not available in an operational 
environment. The CT/AT components are then calculated each 
24 h using the best-track positions relative to the CLIPER 
forecast track. By normalizing the forecast track with 
CLIPPER. the contribution due to climatology and persistence 
1S removed from the forecast. 

There are advantages in using a CT/AT component method 
to measure track forecast errors. The CT component can be 
thought of as turning motion relative to a CLIPER forecast. 
À positive (negative) CT component is defined as a 
best-track position to the right (left) of the  CLIPER 
position at the Same forecast time. Thus, the CT component 
provides a directionality aspect that is not available in 
the forecast error. Similarly, the AT component can be 
thought of as an acceleration/deceleration with respect to 
CETRERS although a change in direction will produce a 
negative AT component even when the translation speed 
remains constant (Neumann and Pelissier, 1981). Ite 
important to note that the CT/AT components are calculated 
relative to a CLIPER reference at a particular forecast 
time. That is, the CLIPER reference is different at 24, 48 
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andes2 h. The importance of this point will become evident 


later. 


D: PFPAST-MOTION CAPEGORIES 

Elsberry and Peak (1986) ftouudl from discriminant 
analysis that the direction and speed of the past  12-h 
movement were significant predictors of future storm motion. 
They divided their data set into five ""past-motion' 
categories based on mash track direction and speed 
(Pie. 2:3)4 These five categories are used in this study. 
The basic assumption of the 5-category specification is that 
the differences in the prior 12-h motion reflect different 
synoptic forcings. For example, there should be a different 


oT sonentdietercing between slow- and fast-moving storms. 


2l 


Likewise, 


storms moving toward the northeast during the past 
12h should experience a synoptic forcing different from 


that of a storm moving toward the southwest or south. 


Cat. 2 


>8 m 








Cat ud (ape 








Northwestward Northeastward 









Cat. 3 


Southward 





Pp Schematic of the five ast-moti t 

based on the past 12-h 5 E e  past-motion categories 
ast pee Categories l and 2 and 

direction GN ENS 3-5) -o£ o cyclone. 


The entire sample of 1357 cases is sorted according to 


the past 12-h motion categories defined in Fig. 20% 


Categories 1 (slow) and 2 (fast) are the speed categories, 


where speed is defined as the translation speed from -12 h 
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to 00 h warning positions. The storms with weak and strong 


steering currents should be isolated in these two 
categories. Categories 3, 4 and 5 are the directional 
categories. For example, storms translating northwestward 
between 2.5 and 8 m/s will be in Category 4. A summary of 


sample sizes in each of the past-motion categories is found 
in Table 1. The largest sample size is found in Category 4. 
The remaining categories have considerably fewer cases. 
Notice that 155 cases are not classified into any of five 
past-motion categories because these cases in the original 


Wilson (1984) were missing a past 12-h warning position. 


Table l1. Sample sizes of the five past-motion categories 
defined by the past 12 h motion of the cycione. 


Category sample Size 
1S2 m/s) 169 
2 (38.0 m/s) 174 
38850-7210? ) 125 
4 (270-340°) 497 
5 (340-80°) 237 
Total r202 


The mean characteristics of the storms within these 
past-motion categories are summarized in Table 2. The mean 
and standard deviation for the speeds of the storms in 
Categories 3, 4 and 5 are approximately the same. Because 


of the large varianée in track direction for the Slow and 


Fast categories, their mean directions are not calculated. 
The mean intensities of cyclones in Categories l, 4 and 5 
appear to be higher. However, the standard deviations in 


these categories are also larger. Therefore, it is unlikely 


that the mean intensities between categories are 
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significantly different. There is also a siight tendency 
for storms in Categories 2 and 5 to be initially farther 
north and west than the other categories, but again the 


standard deviations are large. 


Table 2. Mean (X) and standard deviations (g) within each 
of the five past-motion categories. 


Parameter Cat X O 
Speed (total) l L.6 4.10656 
(m/s 2 WN. 9” 3.5 
i 4.9 1.6 
5.0 main4 
5 ERA 
Speed (x-dir) l -0.5 1.1 
(m/s) 2 -3.2 9.4 
3 =3). 49 93.4 
4 -4.1 1.5 
5 lop eu) 
Speed -dir I 055 e 
E e 
3 de 21.55 
4 2 I.4 
5 4.0 TE 
Directio ? Zu 9 
(dapi 301 18 
5 20 27 
Intensity il 62 26 
(kt) 2 58 24 
3 56.29 
4 69 31 
5 64 26 
Initial 1 13 S 
Longitude 2 136° 12 
(FE) B Pr T5 
4 1h32 aal 
5 136 14 
Thi cao 1 7 6 
Latitude 2 20 8 
N) 3 16 5 
4 17 5 
5 20 7 


The firtst 50 storm tracks within  eaenm of the eve 


categories are illustrated TM: A ION Composite 
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wind fields for each of the categories will be presented 
later. Mast ot the storms in the fast category (Category 2) 
initially moved either toward the northeast or the 
northwest. Therefore, the synoptic flow patterns in this 
fast category might be smoothed when adding wind fields 
associated with the northeast- and northwest- moving storms. 
Notice the general rotation of the past 12-h tracks for the 
storms in Categories 3, 4 and 5 (Fig. 2.5). However, there 
are a variety of future tracks within each category. 

Many of these forecast situations will not have a 
verifiable 72-h position. Within this sample, only about 
255 of the storms have a +72-h position (Table 3). More 
Bao» 30% of the Storms in the fast (Category 2) and 
northeast (Category 5) categories do not even have a +24-h 
position. These differences in storm lifetimes need to be 
considered when interpreting the composites in the following 
section. 

Statistical properties of the CT/AT components within 
each past-motion category are summarized in Table 4. The 
largest sample size is in Category 4 (northwestward) and the 
smallest is in Category 3 (southward). In general, the 
overall means and standard deviations of CT components for 
each time period are smaller than those found by Elsberry 
and Peak (1986), who calculated the CT/AT components 
relative to an extrapolation track. As shown in Table 4, 
the mean CT components are near zero and do not increase 
markedly with time when the CLIPER track is used as a 
reference. The CLIPER normalization also reduces the means 
and standard deviations for the AT components within each 
category and time period compared to those in the Elsberry 
and Peak (1986) study. Therefore, normalization with 
respect to the CLIPER track maybe superior to the use of a 
extrapolation track. Notice that the largest standard 


deviations occur in Category 2 (fast), which is consistent 


23 


Y (KM) 


500 


1500 2000 


1000 


O 





=500 


—1500 —1000 —500 O 500 1000 1500 





O 
o. 
© 
N 
O 
©- 
e 
© 
O 
= 
E 
po 
UM 
> 
oO 
O 
O 
ve) 
o 
O 
Q 
in ' : ! ' 0 
-2000 1500 zIgoO 3900 O 200 1000 1500 
X (KM) 
e. 50 storm tracks from E and 
e e 2n Eee Points plotted are -1 La er eade 
VIL. and 7/2 best-tra@eck NN as 
available. The large dots indicate the -12 h positions. 


26 


Y (KM) 
500 1000 1500 


O 





—500 





-2500 -2000 -1500 -1000 -500 0 500 1000 
X (KM) 


Y (KM) 
500 1000 1500 2000 


O 





—500 


-2500 -2000  -1500 #-1000 —900 Q 509 1060 
E 


1000 1500 2000 


Y (KM) 


500 


~500 





-1500  -1000 -500 0 500 1000 1200 2060 


X (KM) 


Rum. 2.5 ASIA EL. E, except for Category 3 (top), 
Category 4 (middle) and Category 5 (bottom) Broce SEQ 


Dy 
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Tabl : om in 
PEE 2.4 and 2.5 that have verified positions at 24, 48 and 
Time (h) 
Category 24 48 n2 


l («2.5 m/s) 42 33 30 
2 (58.0 m/s) 32 23 16 
3 (80-270?) 39 35 29 
4 (270-3409) 45 39 34 
5 (340-809) 34 22 We 


with the wide variety of track positions at 72 h (Fig. 2.4). 
Some of these storms continued to move rapidly, while others 
Slowed down and began to recurve. 

Histograms of CT and AT provide a graphic representation 
of the data presented in Table 4. For example, the Category 
4 frequency distributions are shown in Fig. 2.6. It will be 
Shown later that it is useful to subdivide the CT/AT 
components into three equal subcategories (terciles) within 
a category. The CT component terciles are to the left, 
center and right of the BIO ERES CI For example, one 
third of the Category 4 CT72 sample is located in the 
central tercile category bounded by CT72 values of -155 and 
108 fof mi. The storm characteristics within each of these 
tercile subcategories will be studied later. A similar 
division into slow, center and fast terciles relative to 
CLIPER is made for the  AT72 components. Since the tercile 
sample sizes become very small in all categories except 
Category 4, conclusions regarding the other tercile 
categories should be viewed as tentative until larger 
samples can be tested. Therefore, only tercile 
subcategories within Category 4 will be examined in this 


study. 
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Table 4. Statistical puc SS of 24, 48 and 72 h CT and 
AT components (n mi) for the past- motion ee oe defined 
in Hig. 2995 The central on third of the tererle sample 
lies between the cutpoints. (N, Sample Size; x, mean; O9 , 
Standard deviation). Overall summary is for 11 five past- 
motion categories of storms during 19/9-1983. 


Tercile Cutpoints 
CT AT 


(UT AT 
Time Cat N DEG INE Ge Left Right SlowwEast 
24 1 135 -19 86 -20 89 -56 14 -50 14 
2 84 5 119 -80 171 -26 46 -116 -4 
Be 740 | ee ccs eee 12-3 1 40 -30 59 
4 367 -4 90 -22 80 -41 35 -54 6 
5 150 1 102 -46 117 -20 44 -81 -3 
overall 810 -2 97 -31 114 
48 1 112 -21 159 -91 204 -75 A A S; 
2 66 -22 25] -22 292 -85 101 -120 36 
3 64 1 224 45 208 -90 92 =455 118 
4 304 -12 181 -60 171 -88 57 -120 0 
5 90 17 201 -140 201 -29 115 -201 -60 
overall 636 -7 203 -54 215 
72 1 95 -17 249 -198 319 -90 70 -283 -34 
2 54 -20 380 -287 585 -166 184 -293 -82 
3 57 -41 36l 3 236 -169 140 -102 109 
4 235 -15 300 -114 250 -155 108 -184 -7 
5 58 -14 292 -160 265 -109 87 -281 -115 
overall 499 -21 316 -151 331 


The mean characteristics of the storms in the CT/2 and 
AT72 tercile subcategories of Category 4 are shown in Table 
5. No significant differences in the past-motion 
characteristics can be identified to explain the tercile 
subcategories. That is, the mean differences are small and 
the standard deviations are large compared to the mean. 
Thus, the storm characteristics at the observation time and 
= 12 tf do not give a. clear indication of the deviations of 
the 72 h position with respect to a CLIPER track. It will 
be demonstrated below that the composite synoptic wind 
fields do provide a discriminator of future motion. 

When using an extrapolated track based on past 12-h 


motion, Elsberry and Peak (1986) Foumce that a majority of 
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Table 5. S (X 
12-h characteristics of Category 
GT or AT terciles at 72 h. 


Means (X) and standard deviations 
Storms 


NOT 

Parameter Cat o Cat 

Speed FLO Left E 1.4 Slow 
(m/s Center 4.8 1.4 Center 

ight 4.8 de 3 Fast 

Speed Wb sin) Left -4.3 1.6 Slow 
il s Center  -4.0 [WM Center 

ight -4.0 1.4 Fast 

Speed yy diz Left Que Las Slow 
(m/s enter Der uS Center 

Right 2.4 I Fast 

Direction Left 299 19 Slow 
degrees Center 300 17 Center 

Right 301 l6 Fast 

Intensity Left 73 31 Slow 
(kt) Center 81 5 Center 

Right 73 33 Fast 

Initial Left 136 Jak Slow 
Longitude Center 134 11 Center 

(7E) Right SS 10 Fast 

Initial Left 14 4 Slow 
Latitude Center 16 4 Center 

MEN) Right 15 3 Fast 
Table 6. Distribution of -2395 EE 4 
cross-track and along-track (AT72) terciles 


AT72 
Slow Center 
Left 25 22 
972 Center 28 B 
Right 24 26 
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the cases in the AT72 fast tercile category were also in the 
CT72 right category. Such a correlation would be expected 
for recurving storms in the extrapolation coordinate system. 
The AT and CT distributions of cases at 7/2 mm, usine a CLIFER 
coordinate system, are given in Table 6. The fast AT72 
cases are fairly well distributed in all three CT72 and AT72 
terciles. In fact, a bias within the CT72 and AT72 terciles 
cannot be demonstrated. This indicates that the tercile 
subcategories, formed using a  CLIPER coordinate system, 
provides a better normalization of the 72-h storm positions 
than does the extrapolation coordinate system used by 
Elsberry and Peak (1986). 
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IB DATA ANALYSIS 


An objective of this study is to demonstrate the 
usefulness of the wind-based EOFs for defining the synoptic 
flow associated with tropical cyclone motion. Shaffer and 
Elsberry (1982) suggested that the environmental flow around 
the tropical cyclone could be described by a small set of 
EOF coefficients based on geopotential height fields. They 
demonstrated that certain eigenvectors were correlated with 


zonal motion of the tropical cyclone, whereas others were 


related to meridional motion. For example, if the EOF 
coefficient 1 is large and positive, the storm initially 
moves northwestward (their Fig. 6). On the other hand, a 


storm tends to move towards the northeast when the EOF 
coefficient ql is small or negative. Peak and Elsberry 
(1986) applied a regression analysis to predict the future 
motion of the cyclone in cross-track and along-track tercile 
categories using EOF coefficients derived from geopotential 
fields as predictors. A second regression scheme was 
developed to predict CT and AT components within past-motion 
categories defined in Fig. 2.3. These sets of regression 
equations provided better track forecasts than those from 
the first set which did not use the past-motion 
categorization. Based on these studies, it seems reasonable 
to examine the synoptic effects in terms of wind-based EOF 
coefficients for cases within categories defined by the past 
12-h motion of the storm. 

To demonstrate that the eigenvectors represent a 
synoptic effect on cyclone motion, composite wind fields are 
calculated using the EOF coefficients. Although compositing 
smooths features particular to individual cyclones, the 
features common to all the cyclones in a synoptic category 


should be evident. By comparing composite fields for each 


2/9 


of the five synoptic categories, it should be possible to 
demonstrate that a small set of the EOF eigenvectors can be 
used to represent differences in synoptic forcing between 
categories. 

The cases in the data set are identified by the 
date-time group and category of past storm motion, so that 
the selection of all the cases in a particular storm motion 
category is easily accomplished. Each case in the sample 
also has a set of (the first) 35 u and v EOF coefficients 
for each of the three pressure levels (i.e., 700, 400 and 
250 mb). For example, the 700 mb EOF coefficients 1 through 
35 may be used to reconstruct the u-component wind field 
(Wilson, 1984) for a given caSe. 

Since the wind-based EOF coefficients (modes) are 
linear, two composite methods are possible. Summation of 
the magnitude of individual EOF coefficients for the cases 
in a category, then divided by the number of cases, would 
yield the mean EOF coefficients. The composite field could 
then be reconstructed by multiplying these mean coefficients 
by the eigenvectors. The alternative method would be to 
reconstruct the fields for each of the cases in a category, 
sum and divide by the number of cases. Both methods will 
yield the same result. The latter method is chosen because 
qnis is the traditional method of compositing for 
meteorological fields. This method will be discussed in 
more detail. 

The composite technique follows the EOF reconstruction 
procedures outlined in Chapter 2. In particular, the wind 
component field is reconstructed gridpoint by gridpoint. 
The normalized wind component for a particular case (j) at 
the gridpoints z(i,j) can be reconstructed from the EOF 
eigenvectors by 

m 
A ` SI i= ee, EE (3.1) 
k=1 
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were we(k,jje ares the EOF coefficients for. the jth wind 
component field and e(i,k) are the appropriate eigenvectors. 
Conversion to a dimensional wind component is accomplished 


by 
AN = 22, ).) Siem bo 121,299... 927, (92 ) 


where b(i) and s(i) are the mean and standard deviation of 
the elements inorow 1 of matrix A. As will be shown later, 
1t is possible to make composite fields using any number of 
Eucrents (e.g. 5, 10, 20, 35, ..., 527). All of the 
variance in the original field would be described by 
including all 527 coefficients. The procedure outlined 
above is repeated for each case in the data set. Composite 
fields are then formed by adding each field and dividing by 


the number of cases in the composite. 


A. PAST-MOTION CATEGORY COMPOSITES 

The composite 700 mb u and v fields for Categories ^ and 
Date shown in Figs. 3.1 and "3.2. The mean u field (Fig. 
3.1) for the northwest-moving (Category 4) storms has a 
Strong easterly component throughout the center of the grid, 
whereas the northeast movers (Category 5) have weaker 
easterly components. The most striking feature in the v 
field (Fig. 3.2) is that the southerly component east of the 
vortex is stronger in the Category 5 composite even though 
the mean intensities of the storms are approximately the 
same in these categories (see Table 2). Tropical cyclones 
moving with a northward bias (Category 5) should have a 
stronger southerly component and a weaker easterly component 
than those moving with a westward bias (Category 4). 
Therefore, the observations in the mean u and v composites 
are physically consistent with the mean directions of 
Category 4 (301°) and Category 5 (20°). Thus, the 35 
wind-based EOFs represent differing synoptic environments 


of the tropical cyclone. 
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The composite véctor wind fields (Fig. 3.3) also provide 
insight into the synoptic differences between the 
past-motion categories. The subtropical ridge is east of 
the vortex in Category 5, but farther to the northeast for 
Category 4. Both composites show a convergence area to the 
southeast of the vortex. The composites tend to smooth any 
midlatitude trough and ridge patterns in the northern part 
of the grid, along with other circulation features that are 
not found consistently in the same location relative to the 
tropical cyclone. However, the mean wind field for Category 
5 gives a stronger indication of a trough to the north of 
the vortex. The cyclonic circulation around the vortex also 


extends over a larger area inthe Category 5 composite. 


Thus, the radial extent of the average vortex in Category 5 
may be greater than that found in Category 4. A similar 
result appears in Chan (1985), who found that the cyclones 


in his northward motion stratification were on the average 
larger than those with a westward motion. 

Subtraction of the composite wind fields from the two 
categories indicates regions with differences in synoptic 
forcing (Fig. 3.4). For example, the largest differences (5 
to 5.5 m/s) are found to the northeast of the vortex. E 
the wind distribution at each gridpoint is assumed to be 
normal, it is possible to determine statistically if the u 
or v component is significantly different at a gridpoint 
using a Student's t-test. The null hypothesis is that there 
is no significant difference in the mean wind components 
between the two categories. If the mean and standard 
deviation at a particular gridpoint yields a t-value of 
greater than 1.96 (95% confidence limit for the sample sizes 
in Table 1), the wind component at that gridpoint would be 
statistically different between the two samples. At least 
one of the two components is significantly different between 
Categories 4 and 5 at 433 of the 527 gridpoints in Fig. 3.4. 
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zero m/s. Dot indicates vortex location. 
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Fig. 3.4. The 700 mb wind difference field using Category 4 
minus 5 based on 35 EOF modes. A pennant is 5 m/s, a barb 
is l m/s, and a half barbs T O TE EE If either the u or 
the v component difference is significanti yaoa e E 
zero (based on the Student s t-test with 95% confidence), a 
dot is plotted at the origin of the difference vector (a 
total of 433 significant points are plotted). Large dot 
indicates vortex location. 
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This strongly suggests that the synoptic forcing differs 
between the two categories. 

The composite 400 mb u component fields (Fig. B 
reveal additional differences between the categories. In 
Category 5, the westerly components extend southward to the 
latitude of the vortex center, whereas a strong easterly 
component is found north of the vortex at 700mb. This 
indicates that the average cyclone found in this category 
experiences a vertical wind shear between 700 and 400 mb. 
Notice that the mean u component at 400 mb in Category ^ is 
again negative to the north of the vortex. 

The 400 mb vector wind fields for Category 4 (Fig. 3.6) 
show the tropical cyclone as an open wave in the easterlies. 
enan (1985) found a similar result from composites of the 
GBA winds for westward-moving cyclones. He pointed out that 
this could be partially explained by the lack of upper-air 
data in the Vicinity of the tropical cyclone. Another 
reason may be the inadequacy of the grid resolution of the 
GBA to describe the smaller radial extent of the tropical 
cyclone at upper levels. Notice that the Category 5 wind 
field reveals a open wave in the westerlies. These 
contrasting features indicate that the vortex in each 
category is experiencing different synoptic forcing. 

Many of the same features discussed above are also found 
at 250 mb. The zonal composites (not shown) have the vortex 
embedded in a belt of easterlies for Category 4 and 
westerlies for Category 5. As at the lower levels, the 
meridional wind composite (not shown) for Category 5 east of 
the vortex has a stronger southerly component than that 
found in Category 4. The total wind fields (not shown) have 
a closed anticyclone east of the vortex for both categories. 

Composite fields for each of the five past-motion 
categories are made using 5, 10, 20 and 35 EOF modes. 


Category 4 (northwestward) will be used as a standard for 
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composites for Category à (top) and Category 5 (bottom). 
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Die. 3.6. As in Fig. 3.3, except for the 400 mb PEE wind 
composites for Category 4 (top) and Category 5 (bottom) . 
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comparison with the other categories because it contains the 
largest number of cases, and the overall mean movement of 
tropical cyclones in the western North Pacific is very 
Similar to the mean movement in Category 4 (Chan, 1985). 
The past-motion categories 1, 2, 3 and 5 are each subtracted 
from Category 4 and the Student's t-test is performed at 
each gridpoint. Table 7 summarizes the total number of 
geridpoints Tor which the u or the v component is 
significantly different between the two categories being 
compared. The number of significant gridpoints varies from 
96 to 491 out of a possible 527 points. The smallest number 
of significant points is generally found in the 4-1 
category, while the largest number of significant points is 
generally found in the 4-5 category. Using 35 EOF 
coefficients for compositing, the differences in the 7/00 mb 
u component field between Categories 4 and 5 are significant 
at j55"of the"327 eridpointSa Similar results can be found 
in the other difference categories, although the number of 
gridpoints with significant differences is not as large. 
Based on the results of Table 7, a statistical basis exists 
for stating that the wind component fields for Categories l, 
2, 3 and 5 are different from those in Category 4. 

Notice that the number of significant points generally 
decreases as the number of EOF modes in the composite 
increases. For example, the category 4-5 difference u 
composite fields at 700 mb decreases from 285 significant 
points for five modes to 229 significant points far, o 
modes. The Student's t-test value is inversely proportional 
to the sum of the sample variances. Since more of the 
variance in the wind field is associated with 35 modes than 
with only five modes (Wilson, 1984), smaller t-values will 
be found. Thus, the chance of the wind components between 
the two categories being significantly different decreases 


with the inclusion of more EOF modes. 
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Mable 7. Number of gridpoints for which the u or the v 
ponent for Categories l, 2, Sor are Significantly 
parrerent from Gategory 4. The number of EOF modes used in 
the compositing is indicated inside the parentheses. 


Number of EOF d 
Level Field Cat (5) (10) (20) (35) 
700 u 4-1 285 279 236 — 229 
4-2 STO 366 334 3072 
4-3 249 244 263 258 
4-5 427 3.998 364 955 
Vv 4-1 297 2:53 248 Zo 
4-2 358 299 239 226 
4-3 324 305 2/3 242 
4-5 342 2:9 0) 280 WE 
400 u 4-1 220 191 169 168 
4-2 NS EZ DON 297 
4 - 3 Z299 260 291 267 
4-5 455 418 420 405 
v 4-1 S27 278 265 247 
4-2 345 312 249 224 
4-3 341 283 234 213 
4-5 535 280 262 240 
250 u 4-1 126 136 101 96 
4-2 428 381 2951 350 
4-3 305 285 284 285 
4-5 491 466 455 458 
V 4-] 241 202 212 n20 
4-2 356 3215 288 269 
4-3 357 Soul 322 223 
4-5 347 276 284 252 
For all the categories, the mean u and v component 


fields at all three levels (not shown) constructed using 5, 
10, 20 or 35 EOF modes reveal only minor differences. The 
magnitude of the wind vectors inthe five-mode EOF wind 
composites are slightly smaller than those in the 35-mode 
EOF composites. This is because more of the variance in the 
Signal is defined when using 35 EOF modes (Wilson, 1984). 
Another minor difference is in the wind direction for light 
winds (less than 5 m/s). This. does not indicate that 
reconstruction of au field for a particular case can be 


done with only five EOF eigenvectors. Wilson (1984) shows 
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that the first 35 eigenvectors contain signal. However, the 
patterns produced in a component field by the 6-35 
eigenvectors are not consistently in the same position 
relative to the vortex. Therefore, these features are 


smoothed during compositing. 


B. TERCILESPATTERN COMPOSITES 

As mentioned in Chapter 2, each of the five past-motion 
categories can be divided into tercile subcategories based 
on either the AT or CT component with respect to CLIPER at 
each forecast interval (24, 48 and 72 h). By comparing the 
composite wind fields of these tercile subcategories, it is 
possible to evaluate further the suitablility of using 
wind-based EOF analysis to define "synopticity" (i.e., basic 
synoptic patterns affecting cyclone motion). Examples of 
the Category 4 storm tracks in the Left, Center and Right 
subcategories are plotted in Fig. 3.7. Category 4 is chosen 
because it has the largest sample size. The track plots 
show a reduction in the scatter of the /2-h positions 
compared with the overall Category 4 tracks in Fig. 2.5. As 
might be expected, storms in the Left subcategory show a 
strong westward zonal displacement. In contrast, the Right 
subcategory has the majority of recurving cases. 

Composite u and v fields are calculated separately for 
storms within the tercile subcategories defined by the CT/2 
components. Pairs of tercile subcategory fields are 


subtracted and at each grid point a Student's t-test is used 


to determine if the tercile composite fields are 
Significantly different at the 954 confidence level. A 
Summary of the test results is found in Table 8. The 


smallest number of significant points is generally found in 
the Left minus Center category. The 400 and 250 mb fields 
are very similar based on these comparisons. Notice that in 
the Left minus Center category the number of significant 


points decreases with elevation, whereas this number 
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increases with elevation in the Center minus Right category. 
This latter feature would indicate that differences in the 
environmental flow between storms moving along a  CLIPER 
track and those moving to the "right o£ a CLIPER. track are 
more significant at the upper levels than at the lower 
levels. Over half of ehe gridpoints (282 pernts) have 
significantly different 700 mb u components for the Left 
minus Right comparison. Since the largest differences 
appear in the Left minus Right comparison, these composite 
fields are chosen for further study. The composite fields 
will be examined to determine if the synoptic patterns 
differ between the subcategories. That is, can the current 
Synoptic patterns be related to the 72-h cyclone position 


relative to a CLIPER reference? 


Table 8. As in Table 7, except the n of different at 
which the wind components are Significantly different is 
computed between two CT72 terciles of Category 4 storms. 


CT72 Tercile Difference Fields 
Level Field Left-Center Left-Right  Center-Right 


700 u 178 232 123 
V 126 150 44 
400 u 94 259 107 
V 33 120 48 
250 u Te 256 149 
V 36 121 56 


The 700 mb mean u fields for the CT72 Left and Right 
tercile subcategories (Fig. 3.8) appear to be very similar 
to that for the overall Category 4 sample (Fig. ld 
However, the extent and magnitude of the easterlies north of 
the vortex in the Left subcategory are larger than those in 
the Right subcategory. This is consistent with the Left 


subcategory storms having a stronger westward component of 


48 


motion. The mean v component fields have a stronger 
southerly component to the east of the vortex in the Right 
mererle category (Fig. 3.9). This feature appears in the 
composite total wind fields as a strong confluence band to 


the southeast of the vortex in the Right category (Fig. 


EE 10). When the two mean wind fields are subtracted (Fig. 
3.11), areas with significantly different forcing become 
apparent. The largest differences are found to the east of 


the vortex, which indicates that the convergence band and 
the anticyclone to the northeast of the vortex play an 
important role in the motion of the cyclones within these 
two subcategories. The composite wind fields at 400 and 250 
mb (Figs. 3.12 and 3.13) show a more intense anticyclone to 


the northeast of the vortex in the Right subcategory than in 


the Left subcategory. The anticyclone to the west of the 
vortex for the Right subcategory is well defined at 400 mb, 
which suggests a break in the ridge, and hence possible 
recurvature. 


Another demonstration that the EOF analysis can be used 
to define synopticity can be made by studying the composite 
wind fields defined within terciles of the AT72 component. 
The AT72 tercile divisions are defined as Slow, Center and 
Fast with respect to the CLIPER forecast. Fifty storm 
tracks for the AT72 subcategories are shown in Fig. 3.14. 
Notice that the Fast subcategory isolates most of the storm 
tracks that translated the greatest distances in the 48-h to 
72-h interval. The Slow subcategory appears to have a large 
number of storms that continued to move northwestward at 
Slow speeds. However, there are a few storms that moved 
left or right of a northwest track. These large deviations 
arise because the  CLIPER forecast positions at 60 and 72 h 
vary with each individual storm based on persistence and 


initial storm parameters. 
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As before, composite u and v fields are calculated for 
each of the tercile subcategories (Slow, Center and Fast) 
defined by the AT72 components. The u or the v composite 
fields are subtracted and at each gridpoint a Student's 
t-test is used to determine if the tercile composite fields 
are Significantly different (Table 9). The smallest number 
of significant points is found in the Center minus Fast 
Subcategory. The number of significant points decreases 
from 700 to 400 mb in the mean u fields and increases from 
700 to 400 mb in the mean v fields. This indicates that the 
zonal EOF modes at 700 mb should better define the 
differences within these categories and the meridional EOF 
modes at 400 mb should better-define the differences within 
these categories. In general, the number of significant 
points is smaller for the AT72 tercile comparisons than is 
found in the CT72 comparisons (Table 8). This suggests that 
less information content exists in the EOF wind components 
regarding the AT72 track components than for the CT72 track 
components. Because the largest differences in composite 
fields are in the Slow minus Fast comparisons, these two 


Subcategories are chosen for further study. 


Table 9. As in Table 7, except the number of dis 
which the wind components are significantly if 
for the AT72 tercile difference comparisons. 


dpoints at 
ferent are 


AT72 Tercile Difference Fields 
Level Field Slow-Center Slow-Fast  Center-Fast 


700 u 140 137 55 
V 94 112 6 
400 u 77 87 14 
V 115 2 49 
250 u 169 169 35 
V 107 123 65 
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The composite zonal fields at 700 mb have a stronger 
meridional gradient in the vortex region for the Fast 
subcategory (Fig. 3.15) than for the Slow subeacecor,. The 
zonal gradient of the meridional wind composite fields 
across the vortex is also larger inthe Fast subcategory 
(Fie. Smee). Notice the stronger southerly component east 
of the vortex in the Fast subcategory. These trends 
indicate that the total winds are stronger within the vortex 
region in the Fast subcategory. These synoptic features at 
700 mb are well defined in the composite wind fields (Fig. 
321798 The extent of the cyclonic circulation is larger in 
the Fast subcategory. Notice also that the anticyclone to 
the northeast is more organized and much closer to the 
vortex in the Fast subcategory. When the wind fields at 700 
mb are subtracted (not shown), the largest differences (3 
m/s) are found west of the vortex. 

Even larger synoptic differences are found at the upper 
levels. The total 400 mb wind composites (Fig. 3. be) 
indicate a more pronounced wave in the easterlies near the 
vortex for the Fast subcategory. The wind difference plots 
at 400 and 250 mb are shown in Fig. 3.19. As expected from 
Table 9, the meridional component accounts for the majority 
of the significant points plotted at 400 mb. That is, those 
wind barbs with significant differences show a more 
north-south orientation: However, the zonal and meridional 
component differences at 250 mb are important in different 
areas. The meridional component tends to be more 
significant in the northern portion of the grid, while the 
zonal component is significant southwest and east of the 
vortex. 

In summary, the composite analysis has demonstrated that 
the EOF modes can be used to describe statistically 
significant differences in the mean wind fields. The 


Synoptic forcing is clearly different between Categories 4 
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and 5, and this forcing is consistent with the mean motion 
of these categories. The same synoptic features are evident 
in the mean wind composites using five and 35 EOF modes. 
Storms in the Left and Right CT72 terciles within Category 4 
are also associated with different synoptic patterns. 
Although the composite patterns in the terciles resemble 
those derived for Category 4, there are noticeable 
differences. Similarly, storms in the Slow and Fast AT/2 
terciles within Category 4 are found to be embedded in 
different synoptic flow patterns. Therefore, the truncated 
EOF coefficients contain synoptic information necessary to 


explain a large part of the future motion of the cyclones. 
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PVTSBESPDEDROCESSTNGOIHE. ONYCM FORBEAST TRACKS 


A. MOTIVATION 

In prior studies on the performance characteristics of 
the JTWC objective forecast aids, the One-Way Tropical 
Cyclone Model (OTCM) has proven to be one of the best 
objective aids.  Tsui (1984) evaluated the performance of 11 
objective aids for western North Pacific tropical cyclones 
dusing 1979-1982. The OTCM had the lowest mean speed and 
forecast errors, and was second best in track error (see 
G@mapter 1, section C for definitions). Elsberry and Peak 
(1986) evaluated the performance of nine forecast aids in 
terms of the CT/AT components relative to a persistence 
track defined by past 12-h movement. Both the CT and AT 
component samples were divided into tercile subcategories. 
The nine forecast aids were evaluated on how well they 
predicted the observed tercile subcategory at each forecast 
interval. In their study, the OTCM was classified as one of 
the top performing aids used by JTWC. 

The forecast philosophy that has developed at JTWC is 
consistent with the results of these studies. Since the 
OTCM has demonstrated consistent skill in identifying track 
and speed changes, the OTCM forecast is used as the primary 
guidance for determining the future track of a tropical 
cyclene (Sandgathe, 1985). The OTCM forecast position at 72 
h is modified in certain forecast scenarios based on other 
objective techniques and known model weaknesses. Two 
examples of modification to the OTCM include adjusting for a 
slow bias on post-recurvature tropical  cyclones, and 
anticipating a more rapid response by the cyclone to 
magilatmtnde forcing than is indicated by the OTCM during 


recurvature Situations. 
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An objective of this study is to develop and test a 
post-processing scheme for correcting the OTCM forecast 
track positions based on the synoptic forcing represented by 
the wind-based EOFs. The technique, to be discussed later, 
should help the JTWC forecaster recognize when the OTCM 
forecast positions need modification. The emphasis here is 
on the OTCM since it has already proven to be a valuable 


forecast aid. 


B. MODEL DESCRIPTION AND DATA AQUISITION 

The current version of OTCM is basically the model 
described by Hodur and Burk (1978). It is a single-grid, 
three-layer, primitive equation model. Ihe coarse grid 
resolution is 205 km. The lateral boundaries have a one-way 
influence from the FNOC global model (Fiorino, 1985). The 
tropical cyclone is initially defined in the model as a 
vorticity bogus. Since the model is dry, the vortex 
circulation is maintained by a specified analytical heating 
Tunce Lon: Since no frictional or other physical processes 
are included, this model is mainly an advective model. 
Modifications to the 1978 version of the model include a new 
method for locating the model grid relative to the initial 
storm position and a stronger storm bogus. A pre-processing 
technique developed by Shewchuk and Elsberry (1978) has also 
been added. The technique (to be discussed later) adds a 
persistence component of motion in the model forecast track. 

The OTCM forecasts from 1979-1983 were provided by T. 
Tsui of the Naval Environmental Prediction Research 
Facility. The 24, 48 and 72-h OTCM forecast positions are 
available at 00, 06, 12 and 18 GMT. The official forecasts 
from JTWC are also available, although not for every time 
the OTCM is run. 
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C.  OTCM CT AND AT COMPONENTS 

The CT and AT components of the OTCM forecasts are 
calculated relative to a CLIPER forecast track as described 
capte Pige 2.2). A Statistical summary of the CT/AT 
components for the OTCM is shown in Table 10. Notice that 
the mean CT components at 48 and 72 h are negative. That 
is, the average OTCM forecast is to the left of the CLIPER 
forecast positions at these times. The negative mean AT 
components indicate that the average OTCM forecast positions 
are behind those of CLIPER. 


Table 10. Mean (X) and standard deviations (0) of CT and 
T components (km lorems eNMUBCM POrecasts du'tinge"19/9-1983 


N = sample size). The CT and AT components are calculated 
with respect to a CLIPER track. 


ET AT 
Time N X O X O 
24 1136 13 154 -126 189 
ZEN 2195343 -167 543 
72 984 -108 522 -152 596 


The frequency distribution of cases defined by the CT 
and AT components at 72 h is displayed in Fig. 4.1. Notice 
Bhsterhue CT72 Xd NPUEtdon is slightly skewed to the left, 
which indicates that the OTCM has a slight tendency to be 


t of a CDBPER C forecast at 472 h. It is somewhat 
surprising to find AT components as negative as -3000 km. 
This situation could occur when the OTCM predicts 


recurvature and the CLIPER forecast is for a continuation 


toward the north. As before, the sample at each forecast 
interval can be divided into CT and AT tercile 
subcategories. For example, the Left subcategory at 7/2 h 


includes all storms with CT components that are smaller than 
-307 km. 
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OO CORRECTION VARIABLES FOR THE OTCM FORECASTS 
Over the years, it has become apparent that dynamical 
models exhibit systematic errors in the forecast track 


relative to the actual tropical cyclone track (Elsberry and 


Frill, 1980). There are a number of causes for such errors. 
First, the models contain an incomplete representation of 
the physical processes in the tropical cyclone. Second, the 


poor resolution of the inner region of the cyclone can cause 


a systematic phase Speed error due to numerical 
considerations. Third, the introduction of a bogus vortex 
can cause systematic directional errors. An overly intense 


vortex may result in excessive poleward components of 


motion. Finally, an important cause for systematic errors 
may be inadequate observations. When observations are 
absent, the wind analysis  reverts toward the first-guess 
field. In the case of the OTCM, this could introduce a 


zonal wind bias since the meridional components due to 
transient features are absent in the first-guess field, 5% 
of which is the monthly climatology. 

In the past, pre- and post-processing techniques have 
been used to minimize the systematic bias in a dynamical 
model. Shewchuk and  Elsberry (1978) demonstrated a 
pre-processing technique for the OTCM. Adjustment of the 
initial wind fields are used to re-direct the predicted 
Brack toward a known direction of movement. A 
post-processing technique developed by Elsberry and Frill 
(1980) integrated the OTCM backward in time to compare with 
the known past positions of the cyclone. This backward 
integration of the model can be rather costly. Peak and 
Elsberry (1983) proposed a backward extrapolation as a 
cost-effective alternative. In this post-processing 
technique, the backward extrapolated positions are compared 
to warning positions. The basic premise is that errors that 
occur in the backward portion of the track may be used to 


help define the expected errors in the forward positions. 
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A goal of this study is to develop a post-processing 
method that adjusts the OTCM forecast positions to reduce 
the systematic errors of the model. The correction method 
is illustrated in Fig mi2: An OTCM forecast position is 
corrected by subtracting the post-processed cross-track 
(CTP) and along-track (ATP) components. A regression 
analysis will be used to develop equations to predict the 
required CTP and ATP adjustments to the OTCM forecast 
positions at 24, 48 and 72 h. 

The CTP and ATP components are defined as 


CTP = OCT - Bem (4.1) 
and 
ATP = OAT - BAT, (4.2) 


where the best-track and CT/AT OTCM components are defined 
as BCT/BAT and OCT/OAT respectively. Thus, the CTP and ATP 
adjusts the OTCM position to the @jobseryed (Dect trace 
Positron. A statistical summary (ome the —CTP andre 
corrections is given in Table ll. A physical significance 
can be applied to the values found in this table. For 
example, amean negative CTP implies the average forecast 
position for the OTCM is to the left of the verified 
Position: A mean negative ATP for the early forecast times 
indicates the OTCM forecast is generally slow, although the 
OTCM is fast by 72 h (i.e., positive mean). At 72 h, the 
combination of a mean negative CTP (left) and a mean 
positive ATP (ahead) indicates that the OTCM may not be 
forecasting recurvature Situations very well on the average. 

The distributions of the desired CTP and ATP corrections 
at 72 h are illustrated inm Fig. 4.3. Although the 


distributions are rather concentrated Detween -700 and 700 
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Table ri; Mean (X) and standard deviations | c } of the 
ost-processing DL CTP and ATP (km or OTCEM 
o e SE- during 1979-1983 (N = sample size). 


CTE ATP 
Time VN X o XAR G 
24 1083 13- 172 -70 186 
48 860 -35 324 212 ^ 83:65 
72 642 -118 508 Tos) 509 


km, there are a number of cases on the tails of the 


distribution where the OTCM errors (CTP and ATP) become very 


large. Hopefully, the errors in these cases should be 
improved in the statistical correction scheme to be 
developed. 


E. REGRESSION ANALYSIS 


Regression analysis is one of the most widely used 


statistical Coors. Neter and Wasserman (1974) give an 
excellent presentation of the theory and method of 
regression analysis. Briefly, it involves using a linear 


combination of predictors to estimate the value of an 
unknown quantity, the predictand. A primary goal of the 
analysis is to choose a set of predictors that minimizes the 
sum of squares of residuals without overfitting. 
l. Predictors 

A set of predictors is sought to assess five factors 
affecting an OTCM forecast: (1) external (to the cyclone) 
physical forcing; (2) previous cyclone movement; (3) current 
cyclone intensity; (4) date; and (5) initial (warning) 
position. In addition, a set of post-processing predictors 
will be used to assess possible errors via a backward 
extrapolation of the track (Peak and Elsberry, Ore) A 
total of 18/ possible predictors are penerated (labile e 
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Four potential predictors are related to the 
observations of the tropical cyclone at the base time. 
These include the Julian date, the JTWC warning latitude and 
longitude, and the maximum sustained wind speed (predictors 
1 to Gane Gable 12). 

The recent motion of the cyclone is an important 
predictor of most statistical forecast schemes. The OTCM 
forecasts are adjusted for a persistence component via the 
Shewchuk and Elsberry (1978) pre-processing scheme. 
However, nine potential predictors representing past motion 
are still included to evaluate OTCM's persistence component 
of motion (predictobse» MN ME; xm» Tab SOSA These 
past-motion predictors are based on JTWC warning positions 
to simulate operational conditions. 

The EOF coefficients derived in Chapter 2 are used 
to represent quantitatively the environmental wind field. 
. Wilson (1984) developed a promising regression track scheme 
using these EOF coefficients. In the present OTCM 
post-processing study, the predictor set includes only the 
first 25 zonal and meridional modes at each pressure level 
(predictors 14 to 163 in Table. 12). Notice that the EOF 
coefficients are identified as Clwnn, where C identifies the 
predictor as a coefficient, 1 is the level (7 for 700, 4 for 
400 and 2 for 250 mb), w is the wind component (u or v) and 
nn is the EOF mode number. It is anticipated thaę the 
synoptic wind forcing represented by the EOF coefficients 
will be key predictors in this regression scheme. 

This study uses the backward extrapolation method of 
Peak and Elsberry (1983) to define a set of six additional 


predictors (Fig. E The +24 h OTCM position is 
extrapolated backwards to define a -24 h position. This 
position is compared to the -24 h warning position. The 


total, meridional and zonal distances (E24, EY24 and EX24) 


from the backward OTCM positions to the warning positions 
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Table 12. BotEDEPST predictors provided in the regression 


analysis. 
Predictor 
Number Name Description 
1 JULDATE Julian date 
2 LAT Warning position latitude 
3 LON Warning position longitude 
4 SPD Maximum sustained wind speed Mc 
5 VX0012 Ponad AA M movement from to 
- al m/ 
6 VY0012 Meridional cyclone movement from 00 
to -12 h (km/hr 
7 V0012 Total aene movement from OO to 
E m 
8 VX1224 AS C epee yoment from -12 to 
- m 
9 VY1l1224 O ES movement from -12 
O - 
10 V1224 Total ee Le. from- l2 to 
-24 B (km). 
JI VX0024 on ES movement from 00 to 
hr 
12 VY00O24 Megadi onal 4 C gans. movement from OO 
w RT) 
13 V0024 mE cyclone movement from 00 to 
14-38 C7U1-25 700 mb coefficients derived for zonal 
modes l to 
39-63 C7V1-25 O0 mb coefficients derived for 
meridional modes l to 
64-88 C4U1-25 400 mb coefficients derived for zonal 
modes l to 
89-113 C4V1-25 400 mb coefficients derived for 
meridional modes 1 to 
114-138 C2U1-25 250 mb coefficients derived for zonal 
modes l to : 
839-163. (6G7V1-25 50 mb coefficients derived for 
meridional modes l to 
164-166 E12 Disb..Ltrom 12 habackward extrapolation 
EUM 3 = CO once Warning pos 1t10n. 
EX12 Meridional and zonal dist. included. 
167-169 E24 Dist. from 24 h backward extrapolation 
EY24, p to -24 h pa ne BOs vel one 
EX24 eridional and zonal dist. included. 
170-172 D2400 Dist. from hi9ibe pae. tien to 00 h 
DY2400, C ee eridional and 
DX2400 zonal dis included. 
173-175 D4800 Dist. from GÈ MOCH position to 00h 
DY4800, A aE E ORAE eridional and 
DX4800 Tzonal ds included. 
176-178  D4824 Dt IIO ie heaUneeposteron to 24 h 
DY4824, OTCM Pos rCion.MeridTonal and 
DX4824 zonal dist. included. 
179-181 D7200 Dist ME 72 h OTCM position to 00 h 
DY7200, Warning p Bos 1b Lore eridional and 
DX7200 zonal dis included. 
181-184  D7248 Dist from 73 h OTCM position to 48 h 
DY7224, OTCM position. Meridional and 
DX7224 zonal dist. included. 
184-187 D7224 Dist. from 72 h OTCM position to 24 h 
DY7224, OTCM position. Meridional and 
DX7224 zonal dist. included. 
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are used as potential predictors. The same procedure is 
used to define E12, EX12 and EY12. However, since the 12 h 
OTCM position is not available, the backward 12 h OTCM 
position is calculated by taking one half of the distance of 
the 24 h position, and projecting this value backward. 

An additional set of predictors (170 to 187 in Table 
12) is used to characterize the forwanmedmoortionso: the OTe 
forecast track py M m It is hypothesized that these 
forward track predictors may explain a portion of the 
variance in CTP and ATP by distinguishing  beteen different 
forecast track orientations. For example, the magnitude of 
the zonal track distance between 72 and 24 h (DX7224) may 
explain some of the variance in CTP at 72 due to recurving 
storms. 

2. Dependent Sample 

It is very important that the regression scheme 
simulate operational conditions. An OTCM forecast to 
provide guidance at 00 and 12 GMT must be initiated from the 
12-h old prediction fields (Fig. 4.5) because the rawinsonde 
reports are not received at  FNOC until several hours past 
the synoptic times. An O ECM prediction at 06 and 18 GMT 
(actually integrated at about 0730 and 1930) is initiated 
with the 00 and 12 GMT analysis fields. This situation is 
chosen to develop a dependent sample for the regression 
analysis. That is, the GBA wind fields from 00 and 12 GMT 
are matched with the 06 and 18 GMT  OTCM predictions. The 
development of a dependent sample in this way ensures that 
the regression analysis scheme will be consistent with the 
current operational version of the OTCM. The 00 and 12 GMT 
OTCM forecasts can not be used because the 12-h predicted 
wind fields are not available [SO this study. 
Unfortunately, the OTCM was either not requested at 06 and 
18 GMT or the forecasts were not archived from 1979-1981, 
because only four 06 and 18 GMT OTCM forecasts are available 
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Pig 4.4 sehemarre 1llustration of OTCM forecast positions 
E NEA +72 _and backward extrapolation positions: (M12 
and M24). Other identifiers are defined in Table 12. The 
displacements identified by dashed lines, along with the 
orresponding zonal and meridional component displacements 
(Gok shown), are used as potential predictors in the 
regression analysis. 
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from 197 90 50] Therefore, the dependent sample consists 
almost entirely of OTCM forecasts from 19602-1703. APtec S 
OTCM forecasts are matched with the wind-based  EOF cases, 
the sample size was reduced to 161 cases at 72 h. Since the 
Sample size is small, withholding cases for an independent 


sample is not possible. 
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3. Prediction Equations 
A linear stepwise regression analysis is chosen to 


predict the CTP and ATP adjustments to the OTCM forecast 


positions. The UCLA Biomedical computer program  BMDP2R 
(Dixon and Brown, 1979) is usedabor thou sys. The 
stepwise regression analysis ensures that a potential 


predictor is not included in the equalicomit et ees brine, 
correlated with any predictor chosen in earlier steps. A 
minimum F-to-enter value of 4.0 is imposed to ensure that 
the De dict oases Significantly correlated with the 
predictand (Dixon and Brown, 19795. 

Two separate regression schemes are developed. The 


first scheme uses only the EOF coefficients as potential 
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Fable 13. Explained variance (R? for caes cIPande-TP for 
aeee Sion scheme COE in which only EQF coefficients sare 
E ode s predictors and ALL in wich all potential 
predictors are permitted. Separate entries are, given for 
Ael imited to a maximum of 5 or 10 predictors (* indicates 
only nine predictors selected). 


2 
Number of Regression R 
Time N Predictors cheme CTP ATP 
24 267 5 COE 0.24 0.19 
10 COE 0.34 0.29 
5 ALL 02533220533 
10 ALE 0.44 0.44 
48 212 5 COE (057 2280.2 7 
10 €OE 0.36 0.40 
5 ALL Dil SO 
10 ALL *0.44 0.47 
72 161 5 COE D. 3590 29 
Fo COE 0.48 0.44 
5 ALL 0.42 0.33 
10 ALL 05855 505.59 
predictors (referred to as COE in later tables). It was 


demonstrated in Chapter 3 that the wind-based EOF 
" coefficients can be used to represent synoptic forcing of 
tropical “cyclone motion. It is hypothesized that these 
coefficients will also explain a significant portion of the 
variance in the CTP and ATP components. The second set of 
regression equations uses all of the potential predictors 
found in Table 12 (referred to as ALL in later tables). 
Comparison of the predictors selected in the two sets and 
the additional skill produced by the ALL set will indicate 
the relative importance of the EOF coefficients. 

The amount of variance explained by the predictors 


is indicated by the R?, which is defined as 

R? = SSR/SSTO = 1- (SSE/SSTO), (4.3) 
where SSTO is the total sum of squares, SSR is ithe 
regression sum of squares and SSE is the residual sum of 


squares. An R? value of 1.00 would indicate that the 


regression equation explains all the variance in the 


ma 


predictand (1.e.. GIP vor Aber The R? values for both the 
CTP and ATP regression sets are summarized in Table 13. The 
amount of variance explained by the equations ranges from 
197 wem 55^. As expected, the regression set that uses all 
of the potential predictors explains the largest amount of 
variance (Table 13). For example, 55% of the variance in 
CTP at 72 h is explained by ten predictors inthe ALL 
scheme. It is extremely encouraging that the largest amount 
of variance is explained by the regression equations at the 
72-h forecast interval. 

The regression equation for CTP at 72 h (CTP72) is 
shown in Table 14. The CTP72 ajustment is obtained by 
summing the product of the regression coefficients (column 
2) and the values of the associated predictors. The first 
predictor explains 18% of the variance in CTP72 (i.e., R? = 
0.18). The second predictor explains an additional 11%, and 
so on. Although 3% of the variance in CTP72 is explained by 
the lOth predictor, the equation is arbitrarily limited to 
this number of variables. Addition of a large number of 
variables may artificially inflate the explained variance 
and yet produce unstable equations when used with an 
independent sample. 

The regression equations for CTP and ATP developed 


using all potential predictors are shown in Tables 15 and 16 


respectively. Corresponding regression equations for the 
coefficients-only (COE) Scheme are found in the Appendix, 
Tables 20 and 21. Notice that the EY1l12 predictor is 


selected early at all three forecast intervals in Table 15. 
The backward extrapolation variable EY12 is a measure of how 
closely the initial forecast track of the OTCM agrees with a 
meridional component of motion based on persistence. That 
is, future errors in the OTCM forecast track are highly 
correlated with the departure of the OTCM from a persistence 


component of motion early in the forecast period. 
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ESI Repression equation for CTP. (km) at 72 h using 
all potential predictors (see Table 12) in the regression 


analysis. R is the amount of variance explained by the 
predictor. Mean T andes nssrdedeviatrfYon (0) tos the 
predictors are included. 
Intercept Sum. 3 
2 = 
Predictor Coefficient R X O 
C7U7 48.4 0.18 -0.5 4.0 
C2V1 24.7 DT T 4.8 
E 1563 0.06 3r) 105 
C4V3 E Ou 0.04 E 4.5 
C4U20 S 0.03 20 LA 
C/U21 46.7 0.03 0.4 1.9 
C2V21 46.2 OROZ 0.2 2 
C / 02 170 Ow 02 B33 Je. 5 
C7V16 D 0.03 BO 1 2E 
D7224 OS 0.03 839 298 
Total DIS 


The most frequently chosen predictors are the 
wind-based EOF coefficients (Table 17). In the ALL scheme, 
» the EOF coefficients comprise between six and eight of the 
ten predictors picked in the regression analysis. The u 
coefficients at 700 mb are most often picked (15 times), 
while the u coefficients at 250 mb are the least often 
picked (one time). It is interesting that the 400 mb u 
coefficients are selected to correct the relative speed 
(ATP) of the OTCM, whereas the 700 mb u coefficients are 
selected most frequently to adjust the directional (CTP) 
aspect of the track. 


Several potential predictors are not included in any 


of the equations. The only purely persistence predictor 
selected is the -12 to -24 meridional motion (VY1224) in the 
CTP24 regression equation. This does not necessarily 


indicate that the pre-processing aspect of the OTCM accounts 
perfectly for departures from a persistence track. The 
early selection of EY12 predictor, which is also related to 
a departure from persistence, probably explains why the 


other persistence-related predictors are not selected. 
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Table 15. Repression equations for CTP (km) at _ the 24, 48 
sil 72 h forecást idees LN (the AC wSeNeme) lheSomdesern 
which the redictors are selected in the equations 1S 
indicated und the parentheses. 


Forecast Interval (h) 


24 48 M2 
Intercept -3.766 -60312 -11.844 
Predictor 
cau” SUSE 2 POZA 2 48.398 l 
CZVI - ge 13. LOTES 265.715 2 
EA 065392 Ge 0.754 (1 1.202 D 
C4V3 - - -26 08 4 
C4U20 - - -77.641 5 
C702 - 32.112905» 46.682 6 
C2V 2 - - 46.211 7 
(UM - - 1740589 8 
C7V16 - - -30.744 (9 
D7224 - - -0.286 10) 
C7V9 7.400 (5) 12.021 (4 - 
C7V23 - 7A 0) IS - 
C4U13 - O to” 7 - 
C4V19 - 14.063 (8 - 
C4V10 - Ld 32/9 - 
EX24 07332 - - 
C4U6 -12.09954 0 : - 
C 7U9 8.449 (4 : - 
C4V2 4.632 (6 - - 
Cd. mo. 539 (7 - - 
C4V11 5) 2000) - - 
C7U16 9.086 ( 95 = E 
4. Verification 
The potential performance of these regression 
equations is tested on the dependent sample. The 


regression-predicted  CTP md ATP components are used to 
modify the OTCM position (referred to as OTCMP in later 
tables). The forecast error is computed by determining the 
great circle distance from the modified OTCM position to the 
best-track position.. The forecast errors for the modified 
and unmodified OTCM forecasts are summarized in Tables 18 
and T2 The -official (ITC) forecast errors for the 
corresponding cases (when available) are also included in 
the tables. The (unmodified) OTCM forecast errors are 
smaller than the official JTWC forecasts (Tables 18 and 19}. 
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48 and 72 h equations for 
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43 
8.3 
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Total 


AS expected, the mean forecast error increases as the 
forecast interval increases. Notice that the mean and 
standard deviations of the forecast errors at all times are 
smaller for the regression schemes (OTCMP) with either five 
or ten predictors than for the unmodified OTCM. À small 
improvement is found in the means and standard deviations 
when ten versus five predictors are used in the regression 
equations. AS expected, the errors for the EOF coefficient 
only-scheme (Table 18) are larger than the errors for the 
other regression Scheme, which includes a PEER predictors 
(Table 19). It is very encouraging that the regression 
schemes developed to modify the OTCM have reduced the mean 
forecast error at 72 h by over 210 km (Table Ti DE ts 
again emphasized that these results are obtained using the 
dependent sample. Independent testing is required to 


determine if the equations are stable. 


Table 18. Mean (X) and standard deviation (o) . of the 
forecast errors (km) for the regression scheme which uses 
EOF coefficients only a with either 5 or 

reddi ctos Also included are forecast errors for a 
once set of unmodified OTCM forecasts and ' a nearly 
homogeneous set of JTWC forecasts. 
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A desirable feature in any forecast scheme is 


consistency, which is indicated by small standard deviations 
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Table 19. As in Table 18, except for regression scheme 
which uses all potential predictors (OTCMP) with either 3 or 
Inooredrcters. 


Forecast Errors 
Number of 


Time N Predictors Method X O 
24 245 = JTWC 206 127 
i = OTCM 203 L2 
" 5 OTCMP NS 95 
10 OTCMP 143 88 
48 190 = JTWC 443 285 
200 = OTCM 386 202 
A 5 OTCMP 302 177 
10 OTCMP 286 148 
72 151 - JTWC 631 400 
161 - OREA 593 oo 
» 5 OTCMP 452 286 
10 OCHE 383 256 

in the forecast errors. Thus, Peers stlenaticantetmak the 


standard deviations for the ALL regression scheme are about 
100 km smaller than those of the unmodified OTCM (Table 19). 
mais result, along with lower mean forecast errors, 
indicates the regresSion equations are extremely successful 
in modifying the OTCM in this dependent sample. 

A schematic representation of the performance of the 
regression scheme using all potential predictors ies 
illustrated in Fig. 4.6. A perfect regression scheme would 
have resulted in all of the cases lying on the  y-x line. 
The correlation coefficient between the observed and 
predicted CTP values is 0.74. The correlation coefficient 
ter ene “APP distribution is 0.73* When there are large 
errors in the observed CTP and ATP, the regression scheme 
can be very successful in adjusting the OTCM to correct for 
those errors. However, there are also times when the 
modification actually degrades the OTCM forecast. 

In summary, two regression schemes are developed. 
The first set uses only the wind-based EOF coefficients as 


potential predictors and the second set uses all the 
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potential predictors in Table 12: The EOF coefficients are 
frequently selected, although the backward extrapolation 
variable EY12 is often the first predictor selected at all 
forecast intervals. This suggests a high correlation exists 
between the OTCM's departure from persistence and future 
track errors. Testing of the regression equations on the 


dependent sample shows the scheme to be very successful in 


reducing the mean forecast errors and the standard 
deviations. The reduction in the standard deviation 
indicates that the modified OTCM forecasts are more 
consistent than the unmodified ones. Independent sample 


testing is not possible since all available cases have been 


used to develop the regression equations. 


F. a CAUTIONSJ FOR THE USE OF THE RECGCRESSTON MODES 

The performance of the regression model on the dependent 
sample is likely to be superior to that from a new sample of 
cases. Such a prediction bias may occur if the predictors 
for the regression model are uniquely related to the 
dependent data cases. The performance of the regression 
model with independent cases depends upon whether conditions 
for the independent case are similar to those used to 
develop the regression model. 

The forecaster should be aware of when to use or not use 
the regression scheme. For example, the forecaster should 
examine rare event situations, such as multiple cyclones, as 
the regression scheme may not perform very well for such 


events. 
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V. CONCLUSIONS AND SUGGESTED RESEARCH 


Using an EOF representation of the wind fields around a 
tropical fcyclone, Wilson 9ST) developed a successful 
statistical track prediction scheme. This study has further 
demonstrated that these EOF coefficients can be used to 
describe synoptic effects on tropical cyclone motions 
Composite zonal and meridional wind fields have been 
calculated for five tropical cyclone categories defined in 
terms of the past 12-h motion of the cyclone. Different 
synoptic flow patterns are found in the environments of 
vortices that were moving northwestward (Category 4) and 
those moving northeastward (Category 5). Statistically 
significant differences in the u or v components based on a 
Student's t-test are found between wind fields in Category 4 
and those in the other four categories. The number of 
significant points over the horizontal grid ranged from 96 
to 49l out of a possible 527 points. The results of these 
t-tests (Table 7) provide statistical supporti LorSE e 
hypothesis that the wind fields for Category ^ storms are 
different from those for storms in the other four 
categories. 

Sets of regression equations are also developed to 
post-process the One-way Tropical Cyclone Model (OTCM). 
Although the backward-extrapolated predictors introduced by 
Peak and Elsberry (1983) are often selected first, the 
wind-based EOF coefficients are frequently selected. The 
schemes are extremely successful in reducing forecast errors 
based on tests with the dependent sample. For example, the 
mean OTCM forecast error at 72 h is 593 km, whereas the mean 
error for the best regression scheme is only 385 km. The 
regression schemes also reduces the standard deviations of 


the forecast errors by as much as 100 km. This eduction 
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indicates that the regression-modified  OTCM forecasts are 
more consistent than the unmodified OTCM. Therefore, it 
appears that the regression scheme has great potential for 


operational applications. 


A. POTENTIAL FOR USE WITH INDEPENDENT DATA 

It has been demonstrated that the EOF coefficients 
provide a large portion of the error reduction in the OTCM 
rack positions. These EOF coefficients are synoptic 
predictors in that they represent the flow fields 
surrounding the tropical cyclone. For a new case, the EOF 
coefficients would be derived using the eigenvectors from 
the dependent sample cases. It is stressed that about half 
of the wind-based EOF coefficients in this study have been 
derived using Wilson's dependent eigenvector matrix. Wilson 
(1984) demonstrated that the eigenvectors derived from the 
dependent sample should be stable since the dependent sample 
size is large (682 cases). However, testing is required to 
determine that the present sample size is sufficient to 
ensure stability. The reader is referred to Shaffer (1982) 
for a detailed example of the appropriate methodology to 
test for stability. 

The regression schemes developed in this study have not 
been tested with independent data. New predictors 
associated with an independent case may differ significantly 
from the dependent sample predictors. When this occurs, the 
performance characteristics of the model are likely to be 
degraded. Therefore, testing on independent data is 
required before the regression scheme can be implemented 


operationally. 


B. OPERATIONAL IMPLEMENTATION 

The regression model developed in this study applies 
only to the 06 and 18 GMT OTCM predictions. No regression 
model is developed for the 00 and and 12 GMT OTCM forecasts, 


89 


which are initialized using 12-h old analysis fields that 
have been slightly adjusted towards climatology in the 
tropics. These adjusted wind fields are not available for 
deriving the EOF coefficients required in this study. 

There are operational limitations for the regression 
scheme developed in this study. The 00 and 12 GMT GBA must 
be available for the EOF analysis. The same storm selection 
criteria described in Chapter 2 has to be used in an 
operational version. Finally, the tropical cyclone must 
satisfy the conditions necessary for making a CLIPER 
forecast (see Chapter 2, section C). 

The operational implementation of the regression scheme 
would be straight-forward. Four operations are required for 
implementation. First: the 25 required EOF coefficients 
must be calculated for each new case. The GBA zonal and 
meridional wind components at 700, 400 and 250 mb must be 
interpolated onto the equidistant grid using the bilinear 
interpolation method (Wilson, 1984). The EOF coefficients 
would then be derived by multiplying the transpose of the 
truncated eigenvector matrix (25 X 527) and the normalized 
observational matrix (527 X 1) for the new case. No 
Significant error is anticipated from normalization using 
the means and standard deviations from Wilson's dependent 
sample. The second step for implementation would be to 
activate the OTCM and CLIPER models, and calculate the OTCM 
cross-track and along-track components with respect to the 
CLIPERCtPQ3CI The third step would be to calculate the 
remaining predictors included in the regression equations. 
Finally, the appropriate predictors would be substituted 
into the set of regression equations to predict the OTCM 
modification values CTP and ATP at 24, 48 and 72 h. Both 
the modified and unmodified tracks would then be transmitted 
to JTWC. As all the information necessary to apply the 


regression modification to the OTCM is available at FNOC, 
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less than one minute of additional computer time would be 


required to also provide the modified OTCM track to JTWC. 


(SUGGESTED RESEARCH 

Although the regression scheme developed in this study 
produces very good results, it may be further improved if 
other storm-related or synoptic predictors are provided for 
the regression analysis. For example, the size of the 
tropical cyclone (as defined by the radius of 30 kt winds) 
could be used. Another example would be to use EOF 
coefficients derived from predicted wind fields, such as the 
24-h forecast from the global prediction model. These EOF 


coefficient predictors might have great potential for 
improving the regression Scheme at longer forecast 
intervals. If the wind fields used to initialize the 00 and 
12 GMT OTCM forecasts were archived, a Similar procedure as 


in this study could be used to develop regression equations 
for modifying these OTCM tracks. 

Since the post-processing procedure proposed in this 
study has shown great potential in reducing forecast errors 
from the OTCM, the same technique could be used to modify 
other forecast aids provided to JTWC. For example, a 
modified Nested Tropical Cyclone Model track could be easily 
implemented using the procedure described in Chapter 4. 

The decision-tree approach introduced by Peak and 


Elsberry (1985) could also be tested as an adjustment 


technique for the OTCM. The algorithm developed by Breiman 
et al. (1984) provides a decision tree with multiple 
branches. Each branch represents a sub-division of the data 
based on the best available predictor for splitting. The 


result is a tree which classifies similar cases and provides 
key values of the classification. The least accurate 
branches are pruned by repeatedly determining the minimum 
errors in the prediction using subsets of the test sample. 


Each branch ultimately ends in a terminal node. A 


91 


statistical summary Of the numper "OT successful and 
unsuccessful classifications in that node from the dependent 
sample provides a measure of the likely accuracy of the 
classification: The probability measure for each terminal 
node would suggest the confidence that the forecaster should 
have in that guidance. It seems likely that the wind-based 
EOF coefficients would provide useful predictors for 
developing such a decision tree. 

In Wilson (1984) and the present study, the wind field 


is decomposed into zonal and meridional components to do a 


scalar EOF analysis. However, a vector EOF analysis of the 
total winds (Legler, 1983) also may be developed to 
determine the possibility of improving the regression 
forecast scheme. In addition, it would be useful to examine 


the effect of rotation of the eigenvectors as suggested by 
Richman. (clase 
In conclusion, the wind-based EOF coefficients provide 


an efficient and accurate representation of the synoptic 


wind forcing oft tropical Wewelione ‘motion: Given this 
scientific basis, application to regression schemes to 
improve tropical cyclone forecasts seems feasible. Further 


research with the wind-based) EOF fields is warranted. 
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